The BCS Pairing Instability in the Thermodynamic Limit 
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The superconducting pairing instability — as determined by a divergence of the two-particle 
susceptibility — is obtained in the mean field (BCS) approximation in the thermodynamic limit. 
The usual practice is to examine this property for a finite lattice. We illustrate that, while the 
conclusions remain unchanged, the technical features are very different in the thermodynamic limit 
and conform more closely with the usual treatment of phase transitions encountered in, for exam- 
ple, the mean-field paramagnetic-ferromagnetic transition. Furthermore, by going to the extreme 
dilute limit, one can distinguish three dimensions from one and two dimensions, in which a pairing 
instability occurs even for two particles. 
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I. INTRODUCTION 

The Bardeen, Cooper, and Schricffer (BCS) theory of 
superconductivity, following the original literature^ is 
typically first presented in textbooks as a proposed vari- 
ational ground state, whose energy is lower than that 
of the corresponding normal state ££■ This is followed by 
a discussion of the excited states, from which, in weak 
coupling at least, a critical temperature is derived, cor- 
responding to the breakup of pairs. 

An alternate view of the transition was put forward by 
Thouless^ who tracked the BCS instability from above 
(in temperature) by monitoring the two particle propaga- 
tor. By including a specific set of processes (denoted by 
"ladders" in the particle-particle channel of the diagram- 
matic version of this formulation) , one finds an instability 
of the normal phase. This approach is also explained in 
many texts^ and will be briefly summarized below. 

Seeing the superconducting transition as an instability 
of the normal state at some finite temperature highlights 
the potential importance of pairing fluctuations that oc- 
cur in the normal state even before the critical transi- 
tion is reached. The possibility of pairing fluctuations 
has been important in the elucidation of the so-called 
pseudo-gap that occurs in the high-T c cupratc materials. 
One school of thought regards the pseudo-gap in these 
materials as a tell-tale signature of pairing fluctuations 
in the normal stated 

There is now an extensive literature on the presence 
of pairing fluctuations in the normal state and on their 
impact on various normal state properties^— However, 
our purpose here is to revisit the simple so-called Thou- 
less criterion for the BCS instability and to reformulate 
it in the thermodynamic limit. We have always found 
it peculiar that this instability is signalled by the ap- 
pearance of two imaginary roots in the denominator of 
the two-particle propagator, whereas other instabilities 
appear to be accompanied by the appearance of a real 
rootJ^ It turns out that the two-particle pairing insta- 
bility is always, to our knowledge, formulated for a finite 
system; in the thermodynamic limit, as we show below, 
the criterion behaves quite differently and much more like 



other instabilities in condensed matter. 

For simplicity we focus on the simplest model that 
exhibits superconductivity (at least at the mean field 
level), the attractive Hubbard model on a hypercubic lat- 
tice in one, two, and three dimensions*^ While consider- 
able attention has been devoted to single-particle prop- 
erties, since these are often most related to the measured 
properties^ the two-particle properties are the ones that 
are key to understanding single-particle properties^ For 
our purposes, we express the two-particle propagator in 
the non-self-consistent ladder approximation^ 
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where xo is the "noninteracting" pair susceptibility 
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and Gocr(k, iu> m ) = \uo m — (ek — fi)] l is the noninteract- 
ing single-particle propagator. Here, k and q are wave 
vectors, and £k is the single electron dispersion appropri- 
ate to tight-binding with nearest neighbour hopping. The 
Matsubara frequencies are defined as iuj m = TrT(2m — 1) 
for Fermions and iv n = i2irTn for Bosons. j3 = (fc^T) -1 
is the inverse temperature, and iV is the number of lat- 
tice sites. All wavevector summations span the entire 
Brillouin zone, and Matsubara sums go over all integers. 
For the BCS instability we can focus on q = and, in 
fact, v n — 0; nonetheless, we wish to illustrate the insta- 
bility by monitoring 172 as a function of (real) frequency. 
The result for xo is 
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where we have now analytically continued the result 
to the upper half-plane (iu„ — > z), and in particular 
for z = v + iS, with 5 a positive infinitesimal. Here 
f(x) = \j(eP x + 1) is the Fermi-Dirac distribution func- 
tion. Equation ^ is the one displayed in textbooks 
and reviews. We show its real and imaginary parts in 
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FIG. 1. Real (a) and Imaginary (b) parts of the noninter- 
acting pair susceptibility at zero wave vector vs. frequency, 
for three different temperatures. Note that poles occur at the 
energies corresponding to two single electron energies. In ad- 
dition, the minimum at zero frequency diverges as the temper- 
ature goes to zero (not evident in (a) because the divergence 
is logarithmic). The figures were produced in one dimension 
with a finite lattice of length 32 sites. The horizontal line at 
0.5 in (a) denotes the value of 1/\U\ for \U\ = It. As the 
central minimum rises above this horizontal line (with de- 
creasing temperature) an instability is signalled, as discussed 
in the text. 



Fig. 1, for a finite (ID) system with TV = 32 at half- 
filling (/i = 0) and for a number of temperatures. The 
chemical potential can either be held fixed or tuned to 
produce a desired electron density, n. 

The use of a finite lattice necessarily produces a series 
of poles, as is evident in Fig. 1(b), where the imaginary 
part of xo is plotted (a small numerical value of 5 is used 
in place of an infinitesimal). These poles correspond to 
the noninteracting case (\U\ = 0); they arise from roots of 
the denominator in Eq. ([3]) with a small negative imag- 
inary part and hence lie in the lower half plane. This 
fact renders them "innocent," since, when substituted 
into the original time-dependent propagator, e~ lvt , they 
will decay away with time. They correspond to normal 
two-particle excitations and result in a density of states 
that extends across the two-particle continuum. In this 
one dimensional case, these energies extend from twice 
the energy of the bottom of the band (— 4i) to twice the 
energy of the top of the band (At). At high tempera- 
tures the excitation energies occur at energies very close 
to those of the noninteracting case, indicated by the ver- 



tical asymptotes in Fig. 1(a); in this case the energies 
are all slightly shifted by the interaction and occur not 
at the asymptotes but at the intersection of the curve 
with the horizontal black line, representing l/|t/|. These 
intersections correspond to the zeros of the denomina- 
tor in Eq. ([T]) and will be referred to as the poles of the 
two-particle propagator. 

A special case, with a distinctive temperature depen- 
dence, is the pair of zeros near v = 0. As the temperature 
decreases, the zero frequency minimum in Fig. 1(a) in- 
creases in value and eventually crosses l/|t/|. As it does 
so, two real roots (with small negative imaginary part) 
become two pure imaginary roots, one of which is in the 
upper half plane. Substitution of this root into e~ wt re- 
sults in an excitation that blows up with time, indicative 
of an unstable normal phase<2i£^ 

We now outline the situation in the thermodynamic 
limit, which, while giving the same physics, looks quite 
different. For q = the discrete sum in Eq. ^ is con- 
verted to an integral over the single-particle density of 
states, g(e)i^ 
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where ±W/2 is the top (bottom) of the single electron 
band (±2i in ID, ±4t in 2D, and ±6i in 3D), and since 
q = we have omitted it from the argument list for xo- 
This integral requires a principal value part, which can 
be done analytically: 
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The integration on the last line is no longer singular and 
can be computed by quadrature. 

In Figs. 2(a) and 2(b) we show the corresponding 
results for an infinite system (in ID) at temperatures 
T = 1, 0.1, and 0.01 (in units of t; hereafter all ener- 
gies will be quoted in units of t). The real part of Xoiy) 
clearly shows a maximum at zero frequency; elsewhere 
there are no positive divergences as they have been in- 
tegrated to a smooth curve in the principal value sense. 
The negative divergences occur at the band edges and 
are due to the divergent single electron density of states 
at the band edges in one dimension. 

As is apparent from the figure, these "band edge" di- 
vergences are present at all temperatures. In fact, for 
the lowest two temperatures shown, the curves are essen- 
tially the same except for the region near zero frequency, 
where the maximum diverges as T — > 0, indicative of a 
superconducting/chargc-dcnsity-wave instability. In fact 



o 



E 



1.00 
0.50 
0.00 
-0.50 
-1.00 
0.50 
0.25 
0.00 
-0.25 
-0.50 



Bulk Limit — 1-D ;; 
u/t = jA- 




Bulk Limit— 1-D 
u/t = 

-8 -6 -4 -2 

v/t 



FIG. 2. Real (a) and Imaginary (b) parts of the noninteract- 
ing susceptibility at zero wavevector vs. frequency, for three 
different temperatures, for the bulk limit in one dimension. 
Note that part (a) in particular looks very different from the 
finite size counterpart in Fig. (la). In particular, the poles 
corresponding to sums of single electron energies are evident 
only in (b). The developing singularity at zero frequency re- 
mains, as is evident in (a). 



at zero temperature the real part is given analytically by 
the following expression (for /! = 0): 
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where v = i>/(4£). The divergence at zero frequency is 
evident in this expression. At finite temperature an ex- 
act analytical expression is not possible, even for zero 
frequency. However, to a very good approximation one 
can obtain Kexoi^ = 0) = ^j log (1.13^). Notice that 
the argument of the natural logarithm is a factor of 2 
larger than what would have been obtained by simply 
approximating the density of states as a constant at the 
chemical potential (// = in this case). 

For a nonzero attractive interaction, an instability is 
signalled by the maximum crossing the black horizontal 
line positioned at 1/|E/|. This signals the onset of an 
instability in a way that is familiar from studies in mean- 
field ferromagnetism,— for example. It appears as though 
two real roots are emerging; in fact a careful analysis of 
Eq. ((T|), using Xo{ v + iS) ~ &o + i c o v , with ao and Cq 
positive real constants [see Figs. 2(a) and (b)] near v « 
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FIG. 3. Real (a) and Imaginary (b) parts of the noninteract- 
ing susceptibility at zero wavevector vs. frequency, for three 
different temperatures, for the bulk limit in three dimensions. 
In 3 dimensions no singularities occur at ±12t, since the single 
electron density of states starts smoothly from zero. The log- 
arithmic divergence in temperature remains at zero frequency. 



shows that as ao traverses 1/|E/|, the same root with 
negative imaginary part becomes a root with a positive 
imaginary part. As in the finite lattice case, then, the 
two particle propagator becomes unstable in time. 

Figure (2b) shows the spectral function, Bo (q, v) = 
— Imxo(q, v + i8)/n as a function of frequency. Aside 
from the asymmetrization, this quantity provides an im- 
age of the single electron density of states. This remains 
true in any dimension, as can be seen from taking the 
imaginary part of Eq. ([5]) : 
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As the temperature approaches zero, the hyperbolic tan- 
gent function simply changes the sign of the density of 
states at the origin. In Fig. (lb) the delta function struc- 
ture was merely providing an image of the finite system's 
discretized density of states. 

As mentioned, this calculation can be done in any di- 
mension, and Fig. (3a) and (3b) show the real and imag- 
inary parts of the pair propagator in three dimensions, 
with nearest neighbour hopping only, at half filling. Once 
again the essential feature is that the real part diverges at 
zero frequency with decreasing temperature; thus, within 
mean field theory, Eq. ([T]) guarantees a transition when 
the real part crosses l/|t/|. Of special note is the lack of 
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FIG. 4. Real part of the noninteracting susceptibility at zero 
wavevector vs. frequency, for three different temperatures, 
for the bulk limit in two dimensions, with nonzero chemi- 
cal potential (a) above the bottom of the band, and (b) be- 
low the bottom of the band. The chemical potential in (a) 
(fi = — 2.5t) corresponds to some low filling of the band. Note 
that the divergence in the real part occurs again at v = 0, 
while the imaginary part of the susceptibility remains pinned 
to zero at this frequency. In contrast, the chemical potential 
in (b) (n — —5. Of) lies below the bottom of the band; this 
means that the real part of the susceptibility diverges at all 
temperatures, and signals a bound state being formed, even 
by two electrons. In 3D this does not occur for a single pair. 



is peculiar to 2D and ID only and shows that any \U\, 
no matter how small, will lead to an instability in the 
extreme dilute limit. This is not the case in higher di- 
mension. For 3D, there is no band edge divergence, and 
it is well-known that an attractive potential must exceed 
some threshold before it will support a bound state for 
two particles. 

Results in other dimensions away from half-filling are 
similar, albeit with differences reflecting the different 
densities of states, as already noted at half-filling, and 
the critical difference just noted regarding band edge 
divergences in 3D vs. 2D and ID. Furthermore, com- 
pletely symmetric results occur for \i = +2.5 (compared 
to /i = —2.5), due to the particle- hole symmetry in the 
problem. 

In summary we have computed the two particle pair- 
ing susceptibility in the thermodynamic limit, in a va- 
riety of dimensions and for any filling. We have shown 
how the BCS instability comes about with decreasing 
temperature and how the nature of the instability more 
closely resembles the one usually discussed in the context 
of mean-field fcrromagnctism. Technically, at the insta- 
bility temperature, a single pole passes from the lower 
half-plane to the upper-half-plane in complex frequency. 
This is in contrast to the finite lattice result, where two 
real roots become two pure imaginary roots, one of which 
leads to the instability. In either case the change at the 
instability signals a two-particle propagator that diverges 
with increasing time. 

The calculations in the thermodynamic limit enable 
one to see the dependence on dimensionality. Some quan- 
titative differences occur due to the very different single- 
particle densities of states in ID, 2D and 3D. However, in 
the extreme dilute limit, the different physics in 3D vs. 
ID and 2D is highlighted in the thermodynamic limit, 
and leads, in a natural way, to the necessity in 3D for 
Cooper's famous calculation^ 



a negative divergence at the band edges in three dimen- 
sions; this affects the extreme dilute limit, and we will 
comment further below. 

Finally, while half-filling gives rise to symmetric re- 
sults, more often than not some other instability inter- 
venes to suppress the pairing instability in this case. A 
result illustrating the pairing instability at a chemical 
potential away from half-filling is shown in Fig. 4, where 
the real part of the susceptibility is shown in two dimen- 
sions for (a) fj, — — 2.5i and (b) fi = —5.0. Note that the 
divergence in the real part persists in (a) at v — as the 
temperature is lowered. The negative divergences associ- 
ated with the band edge discontinuities (in ID and 2D) 
remain, but at a frequency ±W — 2/x. The discontinuity 
at v = — 2\i (5£ for this case) is only present here because 
of the logarithmic divergence in the 2D density of states 
at the origin and is absent for other band structures that 
lack singularities. In (b) the chemical potential is below 
the bottom of the band; now there is a positive diver- 
gence at v — —W — 2/z (here at 2t). This divergence 
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